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We perform a throughout numerical study of the average sensitivity to initial conditions and 
entropy production for two symplectically coupled standard maps focusing on the control-parameter 
region close to regularity. Although the system is ultimately strongly chaotic (positive Lyapunov 
exponents), it first stays lengthily in weak-chaotic regions (zero Lyapunov exponents). We argue 
that the nonextensive generalization of the classical formalism is an adequate tool in order to get 
nontrivial information about this complex phenomenon. Within this context we analyze the relation 
between the power-law sensitivity to initial conditions and the entropy production. 
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I. INTRODUCTION 

Hamiltonian chaotic dynamics is related to the descrip- 
tion of irregular trajectories in phase space. This char- 
acteristic is closely linked to the instability of the system 
and the entropy growth (see, e.g., Q). Chaotic dynamics 
is in fact associated with positive Lyapunov coefficients, 
which corresponds, through the Pesin equality 0], to a 
positive Kolmogorov-Sinai entropy rate 0. However, 
many physical, biological, economical and other com- 
plex systems exhibit more intricate situations, associated 
to a phase space that reveals complicated patterns and 
anomalous, not strongly chaotic, dynamics [J. In many 
of these cases, the system displays an algebraic sensitiv- 
ity to initial conditions and the use of the Boltzmann- 

Gibbs (BG) entropic functional Sbg = — Y^LiPi^Pi 
for the definition of quantities such as the Kolmogorov- 
Sinai entropy rate only provides trivial information. It 
has been first argued then numerically exhibited 
H 0, II IS E3, and finally analytically shown [H [Q 
that under these conditions the nonextensive generaliza- 
tion ^] of the BG entropic form adequately replaces 
these concepts. This generalization has been successfully 
checked for many cases where long-range interactions, 
long-range microscopic memory, some kind of fractaliza- 
tion of the phase space are inherent. Indeed, in the case 
of one-dimensional dissipative maps, following pioneer- 
ing work |14| on asymptotic algebraic sensitivity to ini- 
tial conditions, it has been numerically |S IS E IS E3 
and analytically proved ^f| that the nonextensive for- 
malism provides a meaningful description of the critical 
points where the Lyapunov coefficient vanishes. More- 
over, through the nonextensive entropy it is possible to 
pro ve a remarkable generalization of the Pesin equality 

On the basis of these results, we present here a 
nonextensive approach to the description of complex 
behaviors associated to Hamiltonian systems that sat- 
isfy the Kolmogorov-Arnold-Moser (KAM) requirements 
(see, e.g., ^H). The phase space consists in this case of 



complicated mixtures of invariant KAM-tori and chaotic 
regions. A chaotic region is in contact with critical KAM- 
tori whose Lyapunov coefficients are zero, and a chaotic 
orbit sticks to those tori repeatedly with a power-law dis- 
tribution of sticking times (see, e.g., 0] and references 
therein). It is known that, under these conditions, one 
main generic characteristic feature of Hamiltonian chaos 
is its nonergodicity, due to the existence of a finite mea- 
sure of the regularity islands area. The set of islands is 
fractal, and thin strips near the islands' boundary, which 
are termed boundary layers, play a crucial role in the 
system dynamics. In addition to these processes, there 
are effects that arise solely because of the high enough di- 
mensionality of the system. One of these effects is Arnold 
diffusion [jj]], that occurs if the number of degrees of 
freedom of the Hamiltonian system is larger than two 
PlEl- In general, for the usual case of many-body short- 
range-interacting Hamiltonian systems, the combination 
of these nonlinear dynamical phenomena typically makes 
the system to leave the generic situation of nonergodicity 
and eventually yields ergodicity. 

Along these lines, anomalous effects for the sensitivity 
to initial conditions and for the entropy production have 
already been investigated within the nonextensive for- 
malism |l8j for a paradigmatic model that exhibits the 
KAM structure, namely the well known standard map 
|l9j . In the present paper we focus on two symplectically 
coupled standard maps. The choice of two coupled maps 
is because this is the lowest possible dimension (d = 4) 
of Hamiltonian maps where Arnold diffusion does take 
place, an effect that has a fundamental impact for the 
macroscopic effects associated with the dynamics of the 
system, since it guarantees an unique (connected) chaotic 
sea (see, e. g., |20l l2] | ). 

In Ref. |22J it was numerically shown that, for low- 
dimensional Hamiltonian (hence symplectic) maps, the 
Kolmogorov-Sinai entropy rate K, coincides with the en- 
tropy K produced per unit time by the dynamical evo- 
lution of a statistical ensemble of copies of the system 
that is initially set far-from-equilibrium, i.e., K = K,. 
More precisely, considering an ensemble of N copies of 
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the map and a coarse graining partition of the phase 
space composed by W nonoverlapping (hyper)cells (with 
nonvanishing d hypervolumes) , at each iteration step t a 
probability distribution {pi}i=\,2,...w IS defined by means 
of the occupation number Ni of each cell, pi = Ni/N 
(YliiPi = l)i an d w e have that 



K = lim lim 

t^oo W— >oc 



lim 



(s BG )(t) 



(1) 



where the average (} is taken considering the dynami- 
cal evolution of different ensembles, all starting far-from- 
cquilibrium, but with different initial conditions. In this 
paper, as far-from-equilibrium initial conditions, we will 
consider an ensemble of N copies of the map all randomly 
distributed inside a single (hyper)cell of the partition. In 
this way the initial entropy is zero, since all but one prob- 
abilities vanish. Averages are then obtained by sampling 
the whole phase space changing the position of the initial 
(hyper)cell. 

On the other hand, one can consider the sensitivity to 
initial conditions 



e(x(o),*x(o),t) 



lim 

|5x(0)||- 



*>n*x(o)ir 



(2) 



that in general depends on the phase space initial position 
x(0) and on the direction in the tangent space <5x(0). If 
the system is chaotic, the exponential sensitivity to initial 
conditions defines a spectrum of d Lyapunov coefficients 
{X^}k=i,2,...d coupled in pairs (d being the phase space 
dimension), where each element of the pair is the opposite 
of the other (symplectic structure). By the Pesin iden tity 
we have that J2\( k )>a^ k " > ~ ^ anc ^ the result in |23 
states then that 



E < A(fe) > 



(3) 



where once again () stands for average over different ini- 
tial data. 

Now, when the largest Lyapunov coefficient vanishes, 
Eq. provides only poor information (0 = , to be 
precise), which is not useful for distinguishing, for exam- 
ple, between weak chaoticity and regularity. In order to 
address this issue, we generalize the previous ap pro ach 
in the sense of the nonextensive generalization |l3| of 
the classical (BG) formalism. Specifically, we define the 
q-generalized entropy production (corresponding to the 
q-generalization of the Kolmogorov-Sinai entropy rate) 
as 



— lim lim 

t^oo W— >oc 



lim 



(4) 



where the nonextensive entropy S qs is defined by 

1 - Ei=i Pi 



S a 



q e - 1 



(q e e K; Si = Sbg) , (5) 



and e stands for entropy. In case of equiprobability, i.e., 
p, = I/VP (Vi), the nonextensive entropy takes the form 
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FIG. 1: Phase portrait of two symplectically coupled standard 
maps, Eq. JgJ, with b = 0.5 and ai = a2 = a. Points represent 
the projection of trajectories on different planes. We trace an 
ensemble of N = 4 4 points, with a uniform initial distribution, 
for < t < 50. (a) a = 0: Integrability; (b) a = 0.2: Weak 
chaos; (c) a = 2: Strong chaos. 



Sq c = hlq 



VP, where \n q x = {x^ q - 1)/(1 - q) (with 
lni x = \sxx) is known in the literature as the q-logarithm 
[2^| . Notice that its inverse, the q-exponential, is given 

by 



exp x 



[l + (l-q)x} (1 / < - 1 -^ (ex Pl a; = e x ). (6) 



In this paper we will show that, in situations where 
the largest Lyapunov coefficient tends to zero, a regime 
emerges (that we will discuss in more details later on) 
where the sensitivity to initial conditions takes the form 
of a g-exponential, namely 



Z(t) = [l + (l-q s )\ g J] 



1/(1-9.) 



exp„(A 9s t), (7) 



with a specific value of the index q s < 1 and of the gener- 
alized Lyapunov coefficient \ qs (s stands for sensitivity). 
Correspondingly, a single value of q e < 1 exists for which 
the generalized entropy (JSJ) displays a regime of linear 
increase with time. Under these conditions, we study 
the relation between q s and q e . We will see that, differ- 
ently to what happens in the case of strong chaos (where 
q e = q s = 1), the two indices q e and q s do not coincide 
in general. We discuss the origin of this fact. 
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II. SENSITIVITY TO INITIAL CONDITIONS 
AND ENTROPY PRODUCTION 

We consider a dynamical system with an evolution law 
given by two symplectically coupled standard maps: 



0=2 



b=0.5 



0i(t+l) = pi(t+l)+6 1 (t)+bp 2 (t+l) (modi), 

= pi(t) + (ai/27r)sin[27r0i(i)] (modi), 

2 (i + l) = p 2 (t+l) + d 2 (t) + b Pl (t + l) (modi), 

p 2 (t + l) = p 2 (t) + (a 2 /2TT)sm[2TTe 2 {t)} (modi), 



(8) 



where 9i, Pi £M can be regarded respectively as an an- 
gle and an angular momentum coordinate, a\,a 2 ,b £ 
R, and t — 0,1,2,... . Since the system is sym- 
plectic V(a,6), it immediately follows that it also is 
conservative, i.e., the Jacobian determinant of the 
map is one (det{5[0i(t + l), P i(t + l),6 2 (t + l),p 2 (t + 
l)]/<9[0i(i),pi(t),02(t),p 2 (>)]} = !)■ For the sake of defi- 
niteness, most part of this paper will deal with the setting 
ai = a 2 = a, so that the analysis will be restricted to the 
case where the system is symmetric with respect to ex- 
change 1^2. If the coupling parameter b vanishes, the 
two standard maps decouple. For a generic value of b, the 
system is integrable for a = 0, while chaoticity rapidly in- 
creases with \a\. For intermediate values of \a\, trajecto- 
ries in phase space define complex structures. The phase 
space is non-uniform and consists of domains of chaotic 
dynamics (stochastic seas, stochastic layers, stochastic 
webs, etc.) and islands of regular quasi-periodic dynam- 
ics [l5|. To illustrate these behaviors, in Fig. da-c) we 
show the projection on various planes of the dynamical 
evolution of an ensemble of N — 4 4 points uniformly dis- 
tributed in phase space at t = 0, for a — 2 (strong chaos), 
a = 0.2 (weak chaos), and a = (integrability) . Remind- 
ing the fact that the shadow of a (fractal-like) sponge on 
a wall has a regular geometry, the presence of structures 
as those that appear for a = 0.2 in the projection of the 
ensemble over coordinates planes is a strong hint that 
even more intricate structures are present in the phase 
space itself. In order to further characterize the strongly 
and the weakly chaotic situations, Figs. El and 03 display, 
for different fixed times, the projection of the evolution of 
an initially out-of-equilibrium ensemble, respectively for 
a = 2 and a = 0.2. In all previous figures we have fixed 
b = 0.5. 

Chaotic regions produce exponential separation of ini- 
tially close trajectories. Before entering in the details 
of our analysis, let us consider how the usual, strongly 
chaotic case is described inside the nonextensive formal- 
ism, both for the sensitivity to initial conditions and the 
entropy production. To extract the mean value of the 
largest Lyapunov coefficient and the mean entropy pro- 
duction rate, we analyze respectively the average values 
(\n q t;)(t) and (S q )(t), over different initial conditions. 
Fig. 01 when q = 1, reproduces the results in |22( for 
al = 3, a 2 = 1 and b = 0.5 (strong chaos). Only for q = 1 
both the average logarithm of the sensitivity to initial 
conditions and the average entropy display a regime with 
a linear increase with time, before a saturation effect due 
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FIG. 2: Fixed time evolution of an initially out-of-equilibrium 
ensemble of N = 5 x 10 3 copies of the map ||HJ for the strongly 
chaotic case a = 2, b = 0.5. (a) t = 0; (b) t = 20; (c) t = 50. 
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FIG. 3: Fixed time evolution of an initially out-of-equilibrium 
ensemble of iV = 5 x 10 3 copies of the map (|HJ for the weakly 
chaotic case a = 0.2, b = 0.5. (a) t = 0; (b) t = 20; (c) t = 50. 
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FIG. 4: Average value of the logarithm of the sensitivity to 
initial conditions (a) and average entropy production (b) for 
the map Eq. JHJ in a strong chaos regime [a\ = 3, 02 = 1 and 
b = 0.5). The sensitivity to initial conditions is calculated 
considering the time evolution of the separation of 10 5 couple 
of initial data randomly distributed in the whole phase space 
with an initial separation of order [|<5x(0)[| ~ 1CP 12 . The 
entropy is calculated dividing the phase space in W — 30 4 
hypercells and analyzing the spreading of an ensemble of N = 
6 x 10 6 data initially distributed inside a single hypercell. The 
calculation was repeated 2000 times, choosing the position of 
the initial cell at random in the whole phase space. 



to the finiteness of the total number W of (hyper) cells. 
If q < 1 (q > 1), the W — > 00 curve is convex (concave) 
for both quantities. Notice that the slope of the linear 
increase of (ln^)(t) is smaller than the slope of (SbgH^), 
as, since the phase space is 4-dimensional, there are two 
positive Lyapunov exponents to be considered in Eq. 

We turn now our attention to the weak chaos. Fig. 
shows that, fixing b = 0.5 and decreasing now a\ = 
<22 = a, the mean value of the largest Lyapunov coef- 
ficient tends to zero as a consequence of the fact that 
dynamics is now more constrained by structures in phase 
space. Particularly, for a < 0.8 an initial phase emerges 
for which (ln£)(t) is not linear (see the inset of Fig. [SJ. 
We refer to this phase as the weakly chaotic regime. As 
a tends to zero, the crossover time r (see its definition in 
Fig. 6) between this initial regime and the one character- 
ized by a linear increase of (ln£)(i) increases, so that the 
weakly chaotic phase becomes more and more important. 
The sensitivity to initial conditions in correspondence 
to this initial regime is power-law. In fact, Fig. [7| ex- 
hibits that, for specific values q s (a) < 1, (ln 9s £)(t) grows 
linearly during this regime and starts later to be con- 
vex when the sensitivity becomes exponential, after the 
crossover. For a approaching zero, the crossover time be- 
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FIG. 5: Average value of the logarithm of the sensitivity to 
initial conditions for the map Eq. JSj with b — 0.5 and various 
values of a. The calculation is the result of the analysis of 10 
couple of initial data randomly distributed in the whole phase 
space with an initial separation of order ||<5x(0)j| ~ 10~ 12 . 
The inset reproduces the same axes of the large figure, for 
small values of the ordinate. The solid lines are guides to the 
eye. 



tween this initial regime and the strongly chaotic regime 
(characterized by a linear increase of (ln£)(i)) diverges 
(see inset of Fig. 0). Note that q s depends on a. For 
integrability we have, as expected, g s (0) = 0. 

The weakly chaotic regime can also be detected by 
an analysis of the entropy production, performed as de- 
scribed in section^] In Fig. |S|we take as an example the 
case b = 0.5 and a — 0.2. For the weakly chaotic regime 
the result is similar to Fig. 01b), but the value for which 
the entropy grows linearly is now q e ~ 0.6 instead of 

q e = 1- 

Fig. [5] synthesizes the behavior of q s and q e for two 
different fixed coupling constants, namely b — 0.5 and 
b = 0.2. At variance with what happens in the strongly 
chaotic regime, in the weakly chaotic one q s and q e are 
not the same. In the next section we discuss the origin 
of this difference. 

Another numerical result is that, while the the general- 
ized Lyapunov coefficient appears to be almost indepen- 
dent from a (A 9s = 0.57 and X qs = 0.50 for b = 0.20 and 
b = 0.50 respectively), the generalized Kolmogorov-Sinai 
entropy rate K qc exhibits, within our numerical preci- 
sion, a slight decrease as a increases. We do not have yet 
a clear explanation for this observation. 
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FIG. 6: Time dependence of the time derivative of the average 
g-logarithmic sensitivity for a = 0.05 and b = 0.5 for typical 



values of q. We verify that q s 
arrow indicates the value of r 



0.12 and A„ 
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FIG. 7: Average value of the g s -logarithm of the sensitivity to 
initial conditions for the map Eq. (JHJ with b = 0.5 and various 
values of a. The calculation is the result of the analysis of 10 s 
couple of initial data randomly distributed in the whole phase 
space with an initial separation of order ||<5x(0)j| ~ 10~ 12 . 
q B changes with a (see also Fig. |5J and it was determined 
by a constant derivative condition during the weakly chaotic 
regime. The dotted lines are guides to the eye. In the inset, 
inverse of the crossover time vs a 0,56 . 



FIG. 8: Average entropy production for the map Eq. © 
in a weak chaos regime b = 0.5 and a — 0.2. Phase space 
was divided into W = 30 4 hypercells and the spreading of 
an ensemble of N — 6 x 10 6 initial data was analyzed. The 
calculation was repeated 3000 times, choosing the position of 
the initial cell at random in the whole phase space. We verify 
that q e ~ 0.6 and K qe ~ 1.2; q e changes with a (see also Fig. 
[SJ and it was determined by a constant derivative condition 
during the weakly chaotic regime. The straight line is the 
result of a linear interpolation. 
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FIG. 9: Behavior of q s and q e as a function of a, for b — 0.2 
e b — 0.5. The two parameters coincide only in the case of 
strong chaos (\a\ » 1), when q s — q e = 1. 
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III. DISCUSSION 



To understand the origin of the discrepancy, in the 
weakly chaotic regime, between q s and q e some simple 
geometric considerations are in order. 

Firstly, we notice the following property of the q- 
logarithmic function: 



where a £ 



In, W a = a ln 9 , W , 
and q, q' are related by 



1 



1-q' 



(9) 



(10) 



Notice that, for all a > 1, q' = 1 implies q = 1, whereas 
q' < 1 implies q' < q < 1 . 

Secondly, for simplicity, we consider the case of a 2- 
dimensional phase space where the dynamical evolution 
is symmetric with respect to the exchange of the coordi- 
nates. Let us call = PF (1) <g> VF (1) the total number 
of bidimensional cells composed by a partition of in- 
tervals in each coordinate. We suppose to set a far-from- 
equilibrium initial ensemble inside a single 2-dimensional 
cell and that the dynamical evolution of each trajectory 
in the ensemble implies a spreading of the ensemble in 
phase space. Let W^(t) and W' 2 )(i) be respectively 
the number of occupied 1-dimensional intervals in one 
of the two coordinate axes and the number of occupied 
2-dimensional cells, at time t (^W(O) = W^(0) = 1). 
The relation between (t) and (t) depends on the 
details of the dynamical evolution. Two limiting cases are 

W {2 \t) ocW {1) (t) and W {2) {t) oc [W {1) {t)f . (11) 

The former is realized for instance when there is a pre- 
dominant direction along which the ensemble stretches, 
so that the dynamical evolution of the ensemble produces 
filaments in the 2-dimensional phase space and it is es- 
sentially unidimcnsional. The latter happens for exam- 
ple when the dynamical evolution in the two coordinates 
is decoupled, so that W^(t) is simply the Cartesian 
product W^\t) <g> W^(t) . Now, if we suppose that 
ln 9s W^{t) oc t, for a specific value of the parameter q s 
(we are using the notation q s consistently with the fact 
that £ represents the growth of a unidimensional arc) , by 
means of Eqs. © and (TSJ we have that ln 9e W^ 2 \t) cx t 
for the values 



and q e = 1 — 



1 



(12) 



respectively, for the two previous limiting cases (in Eq. 
(11)). Notice that q s = 1 => q e = 1 for both cases. Of 
course, Eq. (12) can be straightforwardly generalized for 
arbitrary dimensions of the phase space. 

This geometrical analysis applies to the problem of the 
entropy production, strictly speaking, only if the ensem- 
ble evolves according to a uniform distribution in phase 
space (equiprobability) , so that the value of the general- 
ized entropy is given by In, W. Nonetheless, we will see 
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FIG. 10: Behavior of q 3 and q e as a function of b, for the 
integrable case a — 0. The two indices coincide only for the 
special condition b = ±1 (see text for details), when q s — q e — 
0. For almost all other values of b we verify 1 — q e = ( 1 — q s )/2 
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FIG. 11: Fixed time evolution of an initially out-of- 
equilibrium ensemble of N = 5 x 10 3 copies of the map @ 
for the integrable case a — 0, b = 0.5 (first row) and b — 0.2 
(second row). We display the projection of the ensemble on 
the plane (Ox, 62)- 



now that the same analysis is useful to understand the 
dynamical behavior of the model that we are studying. 
We start by considering the effect of the coupling term b 
in the special case a = 0, that corresponds to integrabil- 
ity. Fig. ^1] represents the behavior of q s (b) and q e (b) for 
this case. As expected, q s is zero, thus displaying that 
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FIG. 12: Behavior of a as a function of a, for b = 0.2 (open 
circles) and b — 0.5 (full circles). For a large part of the 
weakly chaotic region a is almost constant and equal to 2. 
The line is a guide to the eye. We do not show the data for 
a > 0.5 because the time interval within which we define q 3 
and q e is not large enough to provide a confident value of 
these quantities. 



the sensitivity to initial conditions increases linearly. On 
the other hand, the behavior of q e is understood if we 
analyze in some details the dynamics. We have that pi 
and p 2 are just conserved along any trajectory and all the 
interesting dynamical evolution take place in the (0\, 6 2 )- 
plane by means of the iteration laws 



h(t + i) = pi(o) + e 1 {t) + b P2 (o), 

h{t+l) = P2(O)+0 2 (t)+b Pl (O). 



(13) 



As it is apparent, the two coordinates are decoupled and 
for b ^ ±1 the growth of the ensemble is 2-dimcnsional 
as confirmed by Fig. El (first row). Under these condi- 
tions, we have W (i) (t) oc W [2) {t) oc [W (1) (t)} 2 , so that, 
applying the previous analysis, we obtain q e = 0.5. For 
b = ±1 the maps further degenerates, since now we have 
one more conserved quantity, i.e., 



9i(t+l)T0 2 (t + l) = 9i(t)T0 2 (t). 



(14) 



This implies W^(t) oc W^(t) oc W^(t) (see Fig. EI 
second row), and we obtain q e = q s = 0. Around b = ±1 
there is then a transition from q e = 0.5 to q e = 0, as it is 
exhibited in Fig. ITU1 

We can consider now the generic weakly chaotic case 
(0 < \a\ << 1). In Fig. El w e estimate the (below de- 
fined) factor a that connects the 4-dimensional analysis 
performed by the entropy to the one-dimensional inspec- 
tion provided by the sensitivity to initial conditions, by 



means of the relation (see Eq. I|10|) ) 
_ 1 - Qs 



1 



(15) 



For both b = 0.5 and b — 0.2, small positive values of 
a yield a ~ 2 , as it is for the integrable case a = 0. 
This exhibits an essentially 2-dimensional growth of the 
portion of the phase space occupied by the ensemble. 

At this point we can advance a conjecture for the more 
general case of a phase space that is not symmetric under 
the interchange of coordinates. This conjecture will allow 
for a better understanding of the reason of the discrep- 
ancy between q s and q e . In the case of strong chaos in a 
symplectic nonlinear system, the growth of the hypervol- 
ume containing the initially out-of-equilibrium ensemble 
is essentially d/2-dimensional (d being the phase space 
dimension of the system). In fact, if {\^}k=i,2,...,d/2 
with \W > VA; is the spectrum of positive Lyapunov 
coefficients, and we call £w the sensitivity to initial con- 
ditions associated to the direction defining \w , we have 

£ (1) (*K (2) W ••• e (d/2) 

oc exp(A (1) i) exp(A (2) i) 

/ d/2 

oc exp 



W (d \t) 



(16) 

'(*) (17) 
exp(A (rf/2) iI18) 

(19) 



We remind that W^(t) and W^ d / 2 \t) in Eq. (16) re- 
spectively are the numbers of occupied o?-dimensional and 
rf/2-dimensional hypercells. Taking the logarithm of rela- 
tion (19) we obtain essentially the Pesin equality, namely 



d/2 



(20) 



k=l 



where \nW {d \t) ~ Kt. 

If we now assume that the sensitivities to initial con- 
ditions along the various directions of the expansion 
are ^^-exponential power-laws (which implies that the 
largest Lyapunov coefficient vanishes) 



(21) 



the previous relations (16-19) transform, for very long 
times, into 

W {d \t) oc W (d/2 \t) (22) 

oc C (2) W ••• t W2) (t) (23) 

OC *V(1-3 (1) ) t l/(l-9 (2 >) ... W(l-9 <d/2) )(24) 



(25) 



and, since ln ?c W^(t) is linear with time for q e (i.e., 
S qe oc [W^(t)] 1 -i- oc t, hence W^(t) oc t q ">) we 
obtain 



1 d/2 1 
1 - Q, f^l-i 



k=l 



(26) 
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Note that, for < 1 Vfc, we have q e > g( fc );Vfe. Note 
also that for gW = q s Vfc Eq. J2SJ| reduces to Eq. (IT5f) 
identifying a = d/2. 

Since Eq. (|26[) relates an information associated with 
the entropy production to an information associated with 
the sensitivity to initial conditions, it plays, for weakly 
chaotic dynamical systems, the role played by the Pesin 
equality for strongly chaotic ones. Working out a gen- 
eralization of traditional methods for the calculation of 
the complete spectrum of the Lyapunov coefficients (e.g., 
[24[), it might be possible to numerically verify the va- 
lidity of Eq. (HJ 

Eqs. (|2"Uj) and JSBJl can be unified as follows (t » 1) 



d/2 

exp 9e (^ 9c t) ^ [exp^*) (A^2, t) 
fe=i 



(27) 



This relation hopefully is, for large classes of dynamical 
systems yet to be qualified, a correct conjecture. If so 
it is, then it certainly constitutes a powerful relation. 
Indeed, let us summarize some of its consequences: 

(i) If the system is strongly chaotic, i.e., if it has pos- 
itive Lyapunov exponents, we have q e — q s = 1 and Eq. 
(20) holds; 

(ii) If the system is weakly chaotic (hence its largest 
Lyapunov exponent vanishes), we have that Eq. (26) 
holds; 

(iii) If the system is weakly chaotic and there is only 
one dimension within which there is mixing, then, not 
only q e = q s , but also 



K„ 



(28) 



as already known for unimodal maps such as the logistic 
one [Hill; 

(iv) If the system is weakly chaotic and there is isotropy 
in the sense that q^ = q s (Vfc = 1, d/2), we have that 



f 



d/2 

l—q e l-q s 
The special case q s = yields 



Qe = 1 -d- 



(29) 



(30) 



It is suggestive to notice that, discussing a Boltz- 
mann transport equation concerning a fluid model with 
Galilean-invariant Navier-Stokes equations in a d/vs- 
dimensional Bravais lattice, Boghosian et al j2|| ob- 
tained q e = 1 — 2/dws, and that, for a lattice Lotka- 
Volterra fi^y-growth model, Tsekouras et al [27| obtained 
q e = 1 — l/d^v (NS and LV stand respectively for 
Navier-Stokes and Lotka- Volterra) . 



IV. CONCLUSIONS 

We have discussed the (average) sensitivity to initial 
conditions and entropy production for a 4-dimensional 



dynamical system composed by two symplectically and 
symmetrically coupled standard maps, focusing on phase 
space configurations characterized by the presence of 
complex (fractal-like) structures. Under these conditions, 
coherently with previous 2-dimensional observations [18| , 
we have detected the emergence of long-lasting regimes 
characterized by a power-law sensitivity to initial condi- 
tions, whose duration diverges when the map parameters 
tend to the values corresponding to integrability. While 
the classical BG formalism characterizes these anoma- 
lous regimes only trivially by means of the Pesin equality 
(we have K = and \W = Vfc in Eq. ©), the 
nonextensive formalism fl3T| .through a generalization of 
the Pesin equality (see also [EEI); provides a meaning- 
ful nontrivial description for these regimes. Specifically, 
we have shown that, during these anomalous regimes 
(here called weakly chaotic) the average sensitivity to 
initial conditions is a ^-exponential power-law ©, with 
q s < 1, and that the corresponding entropy production is 
(asymptotically) finite only for a generalized g e -entropy 
(0, with q e < 1. We have discussed the relation between 
q s and q e , both numerically and analytically and we pro- 
pose an appealing generalization of the Pesin equality, 
namely Eq. (|27(l . By means of the difference between q s 
and q e we obtain useful information about the dimension- 
ality which is associated to the spreading of the initially 
out-of-cquilibrium dynamical ensemble. If this dimen- 
sionality is equal to one, we have the result q e = o, and 
Kq,, — \„ ! a s already conjectured p| and proved [j2 for 
the edge of chaos of unimodal maps. 

The present results concern low-dimensional conserva- 
tive (Hamiltonian) maps, but the scenario we have de- 
scribed here may serve as well for the discussion of anal- 
ogous effects arisin g w hen many maps are symplectically 
coupled (see, e.g., |2j|), and even for the case of many- 
body interacting Hamiltonian systems. Short-range in- 
teractions would typically yield strong chaos, and long- 
range interactions may typically lead to weak chaos (con- 
sistently with the results in [281 29]). In fact, for iso- 
lated many-body long-range-interacting classical Hamil- 
tonians there are vast classes of initial conditions for 
which metastable (or quasi-stationary) states are cur- 
rently observed, later on followed by a crossover to the 
usual Boltzmann-Gibbs thermal equilibrium (see, for in- 
stance, H3)- The duration of the metastable state di- 
verges with the number N of particles in the system, 
in such a way that the limits N — > oo and t — > oo do 
not commute. Since it is known that such systems may 
have vanishing Lyapunov spectrum, it is allowed to sus- 
pect that the scenario in the metastable state is similar 
to the weakly chaotic one described in the present pa- 
per, whereas the crossover to the BG equilibrium corre- 
sponds, in the system focused on here, to the crossover 
to the q e = q s = 1, strongly chaotic regime. These and 
other crucial aspects can in principle be verified, for in- 
stance, on systems with very many (and not only two, as 
here) symplectically coupled standard maps. They would 
provide insighful information about macroscopic systems 
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and their possible universality classes of nonextensivity ledo, as well as partial financial support by Pronex/MCT, 
(characterized by index(es) q). Efforts along this line are Faperj, Capes and CNPq (Brazilian agencies), 
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